Thermal effects on chaotic directed transport 
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We study a chaotic ratchet system under the influence of a thermal environment. By direct 
integration of the Lindblad equation we are able to analyze its behavior for a wide range of couplings 
with the environment, and for different finite temperatures. We observe that the enhancement of 
the classical and quantum currents due to temperature depend strongly on the specific properties 
of the system. This makes difficult to extract universal behaviors. We have also found that there 
is an analogy between the effects of the classical thermal noise and those of the finite h size. These 
results open many possibilities for their testing and implementation in kicked BECs and cold atoms 
experiments. 
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I. INTRODUCTION 



The first theoretical constructions related with di- 
rected transport have been formulated in an early work 
by Feynman [1]. This opened a field whose relevance 
and activity has been increasing since then. The moti- 
vation is twofold, in the first place several fundamental 
questions have originated from these ideas and have been 
only partially answered [2]. As a consequence, the sub- 
ject grew into a major new field of statistical physics. 
On the other hand, the wealth of possible applications 
has determined the emergence of a very broad area of 
research in physics. Ratchets, generically defined as pe- 
riodic systems having a dissipative dynamics associated 
with thermal noise and unbiased perturbations (driving 
them out of equilibrium), can be used to model a wide 
range of different phenomena. In order to give just a 
few examples we can mention molecular motors in bi- 
ology [3], nanodevices (rectifiers, pumps, particle sepa- 
rators, molecular switches and transistors) [4], and cou- 
pled Josephson junctions [5]. On the other hand, there 
is a great interest in the theoretical description and ex- 
perimental implementations of cold atoms subjected to 
time-dependent standing waves of light. They play a cen- 
tral role in many studies of the quantum dynamics of 
nonlinear systems (dynamical localization, decoherence, 
quantum resonance, etc.) [6]. In particular, the so-called 
optical ratchets, i.e. directed transport of laser cooled 
atoms, have been successfully implemented in this sort 
of experiments [7]. 

The appearance of a net current can be classically ex- 
plained by the necessary condition of breaking all spa- 
tiotemporal symmetries leading to momentum inversion 
[8]. This, along with Curie's principle indicates that 
transport should be present, given that it is not for- 
bidden. As an example of this situation we can men- 
tion Hamiltonian systems (with necessarily mixed phase 
spaces) where a chaotic layer should be asymmetric [9]. 
In the more general dissipative case, chaotic attr actors 
having this property are necessary [10]. It is usual that 
the same principle translates almost directly into the 



quantum domain [11], and very similar behaviors arise. 
But sometimes quantum mechanics introduces new ef- 
fects [12] and the results depart from the classical ones. 
For example, the efficiency of a forced thermal quantum 
ratchet has been calculated in [13]. In that work the 
authors find that the quantum current is higher in com- 
parison with the classical one at the lowest values of the 
temperature. As this parameter increases the discrepan- 
cies diminish and finally they become negligibly small. 

Recently, a quantum chaotic dissipative ratchet has 
been introduced in [14]. In this example directed trans- 
port appears for particles under the influence of a pulsed 
asymmetric potential in the presence of a dissipative en- 
vironment at zero temperature. The asymmetry of the 
quantum strange attractor is at the origin of the quan- 
tum current, in close analogy with what happens at the 
classical level. Indeed, this work provides with the pa- 
rameters needed for a possible implementation using cold 
atoms in an optical lattice. For a somewhat similar dy- 
namics, the case of weak coupling and low temperature 
has been studied in [15]. In the present work we extend 
the study of [14] for a wide range of couplings with the 
environment and different temperature values. We have 
verified that there is a strong dependence of the current 
behavior on the coupling strength. If we compare this 
with the results found in [13], for instance, we could not 
find a generic enhancement of the quantum current for 
finite temperatures. Instead of that, we could identify 
a close quantum-to-classical correspondence when con- 
sidering thermal effects only at the classical level. In 
fact, we have found that the finite H effects on the quan- 
tum current are analogous to the influence of the classical 
thermal fluctuations on the classical transport. 

In the following we describe the organization of this 
paper. In Section II we present our model for the sys- 
tem and for the environment, explaining the methods we 
have used to investigate the current behavior. In Sec- 
tion III we show the results where the roles of h, the 
coupling strength and the temperature are analyzed in 
detail. Finally, in Section IV we summarize and point 
out our conclusions. 
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II. THE SYSTEM AND THE ENVIRONMENT 

In this Section we describe the approach used to model 
the system plus the environment. We study the motion 
of a particle in a periodic kicked asymmetric potential 
given by 

+ CO 

V(x,t) = /c[cosx+a/2cos(2x + $)] ^ S(t-mr) (1) 

TYl— — OO 

where r is the kicking period, k is the strength of the 
kick, and a and are parameters that allow to introduce 
a spatial asymmetry [14]. The effects of the environment 
are taken into account by means of a velocity dependent 
damping and thermal fluctuations. At the classical level, 
these ingredients are incorporated in the following map: 

J n = Tn + k[sinx + a sin (2x + <£)] + £ 

yx = x-\-rn ^ ' 

In these expressions, n is the momentum variable con- 
jugated to x and T is the dissipation parameter, with 
< T < 1 . The thermal noise £ is related to T, accord- 
ing to < £ 2 >= 2(1 — r)/c#T, where ks is the Boltzmann 
constant and T is the temperature, making the formula- 
tion consistent with the fluctuation-dissipation relation- 
ship. By performing the change of variables rn — > p, 
rk if, and r£ £ (where < f 2 >= 2(1 - Y)k B f and 
T = r 2 T), we can eliminate the period from the classical 
expressions, and define the new map 

( p = Tp + K[smx + asin(2x + $)] + £ ^ 
\ x = x + p 

In the quantum version of the model the system Hamil- 
tonian is given by 

H s = n 2 /2 + V(x,t) (4) 

where the quantization has been performed in such a way 
[16] that x — > x, n — > n = —i(d/dx) and ft = 1. This 
amounts to saying that, being [f,p] = zr, there is an 
effective Planck constant given by ft e ff = t. Then, the 
classical limit corresponds to ft e ff —> 0, while keeping K = 
h e ffk constant. 

In order to incorporate dissipation and thermalization 
to the quantum map we consider the coupling of the sys- 
tem to a bath of non interacting oscillators in thermal 
equilibrium. The degrees of freedom of the bath are elim- 
inated introducing the usual weak coupling, Markov and 
rotating wave approximations [17]. This leads to a Lind- 
blad equation in action representation 

9 = -i[Hs,p] + 

+ g^J^ 

+ ^^n7^^ n ,T)n7j^ nS T){[Lt,pL n/ ] + [Ltp,L n ,]} 

(5) 



where the frequencies £7 n = n + 1 /2 are the energy differ- 
ences between two neighboring levels of the rotator. The 
population densities of the bath found in Eq. 5 are given 

by 

n th( U >T) ~ exp(hoj}k B T)-\ 

(6) 

n+>,T) = n t ->,T) + l 

The system operators L n = + 1(| n )(n + 1 1 + | — 
n)(—n — 1 |) describe transitions towards the ground 
state of the free rotator. Requiring quantum to clas- 
sical correspondence at short times we fix the coupling 
constant g = — ln(l — V). For T = we recover the mas- 
ter equation used in [14] for the pure dissipative case. 
For finite temperature, the last term in Eq. 5 describes 
the thermal excitation of the rotator through absorption 
of heat bath energy. Finally, Eq. 5 will be integrated 
numerically without further approximations. 



III. RESULTS 

Since we are interested in chaotic transport, through- 
out the following calculations we will use the set of pa- 
rameters given by K = 0.7, = 7r/2, and a = 0.7. In 
the Hamiltonian limit, this case shows no visible stability 
islands in phase space. We have first studied some classi- 
cal aspects of our system at zero temperature, beginning 
with the bifurcation diagrams in terms of p as a function 
of the parameter T G [0,1] (see Fig. 1(a)). It should 
be mentioned that the chaotic attractors set in very fast 
for small T. But this is not the case for larger values 
where the transient times can be very long. For this 
reason, we have calculated these diagrams with the last 
5 x 10 3 iterations of the map, after the first 1.4 x 10 5 have 
been discarded. We have randomly taken 5 x 10 3 initial 
conditions inside the unit cell (x G [0, 7r), p G [— 7r,7r]). 
Several regular and chaotic windows alternate. The for- 
mer are characterized by simple attractors (stable fixed 
points of the dissipative map) and the latter are domi- 
nated by chaotic attracting sets. The width in p of these 
sets grows as dissipation weakens (r — > 1). 

The main quantity characterizing transport is the cur- 
rent J {t) = (pt), where () stands for the average taken on 
the initial conditions, and p t is the moment after the tth 
iteration of the map. Since bifurcations that change the 
shape of strange attractors also play a role in determin- 
ing the values of the asymptotic current we restrict our 
analysis to the chaotic window approximately located at 
T G [0.68,0.78] indicated by an arrow. The inset of Fig. 
1(a) displays the asymptotic current Joo = lim^oo J(t) 
as a function of T (circles stand for T = 0). At this 
T range, 100 kicks are enough to reach the stationary 
behavior, independently of the initial distributions . As 
pointed out in [14], dissipation induces an asymmetry of 
the strange attractor which is responsible for the directed 
transport. On the other hand this dissipation mechanism 



3 



contracts phase space and makes the higher energies inac- 
cessible for the system. The final value of (p) results from 
the interplay between both effects. In fact, increasing dis- 
sipation does not increase the transport and the largest 
values of (p) are obtained for the lowest values of dissipa- 
tion, i.e. T > 0.9 (for example = 5.78 for T = 0.97). 
The minimum current in this window is reached for an 
intermediate value of dissipation (T = 0.74). 

We then consider the case of finite temperature. The 
bifurcation diagram corresponding to T = 0.05 is shown 
in Fig. 1(b). It is clear that the effect of temperature 
consists of smoothing the finer structure of the chaotic 
attractors that is present for the smallest values of T. 
Even for this extremely low value of T the detailed fea- 
tures have almost completely disappeared, with the ex- 
ception of the black lines corresponding to the highest 
values of the density distributions. The other very inter- 
esting effect is that temperature erases the regular win- 
dows allowing for a continuously chaotic behavior. This 
could be of much relevance in obtaining large ratchet cur- 
rents without the need for an extremely fine tuned, weak 
dissipation [18] (this will addressed in future studies [19]). 
As shown in the inset, low temperatures (diamonds cor- 
respond to T = 0.05) lead to a noticeable enhancement 
of the asymptotic current (around 30% for T = 0.74). 

We now turn to compare the classical and quantum 
behaviors. Firstly, we analyze the currents (which in the 
quantum case is given by J(t) = tr(pp)). In Fig. 2 we 
display as a function of T, for three different values 
of T and h e ff. At the classical level low temperatures lead 
to an enhancement of the current for intermediate values 
of the dissipation (see Fig. 2 upper and middle panels 
corresponding to T= 0.7, 0.75) . In the case of weak dis- 
sipation ( T = 0.9 in the lower panel of Fig. 2), which 
displays larger values of , the effect of thermal noise is 
negligible. For higher temperatures, the thermal effects 
reduce (p) in all cases. This can be interpreted as follows: 
thermal noise reduces the energy loss caused by dissipa- 
tion (with no kicks, the system would attain a Boltzmann 
distribution), so higher energies can be reached, in com- 
parison with a pure dissipative process. But since this 
diffusion also tends to blur the asymmetry of the strange 
attractor, the two effects compete and transport has a 
maximum for low values of T, and then decreases. 

At the quantum level we can clearly see that the pre- 
viously described thermal enhancement is not generally 
present, at least for the ft e ff values we have considered. 
It is important to note that these values are consistent 
with experimental implementations using for example, 
cold sodium atoms in a laser field having a wave length 
A = 589nm (for more details see [6, 14]). However, 
we observe a very slight growth of the current for the 
^eff = 0.055 case with T = 0.9 (see the blue triangles in 
the lower panel of Fig. 2)indicating that the tempera- 
ture dependence of the current is very sensitive to the 
particular dynamics of the system. For a different ex- 
ample, an enhancement of the quantum transport has 
been observed [13], hence it is difficult to extract uni- 



p 




Figure 1: Bifurcation diagrams in terms of p, as a function 
of the parameter F. We show the last 5 x 10 3 classical p val- 
ues corresponding to 5 x 10 3 random initial conditions taken 
in the interval (x G [0, 7r), p £ [— 7r,7r]), and after 1.4 x 10 5 
steps of the map. In panel (a) T — (inset shows the asymp- 
totic current Joo as a function of F for the indicated irregular 
window), while in (b) T = 0.05. 



versal behaviors. For T = 0.7 the quantum current is 
lower than the classical one for all temperatures, as al- 
ready pointed out in [14], but only for T = 0. The same 
happens in the case of weak damping T = 0.9 (neverthe- 
less, for ft e ff = 0.055 and T > 0.1 both currents coin- 
cide). The case F — 0.75 shows a different behavior. At 
T = the quantum currents (for any of the ft e ff values 
we have considered) are larger than the classical ones, 
that is, there is an enhancement due to the finite size 
of h e R. Also, the maximal quantum current corresponds 
to h e ff = 0.165, and this is valid for all the tempera- 
tures shown (for T = 0.05 it is still greater than the 
classical one). For larger quantum coarse-graining (see 
^eff = 0.494) quantum currents decrease. In this sense, 
it seems that the effect of quantum fluctuations on the 
quantum directed currents is analogous to the effect of 
those of thermal origin on the classical ones. Small fluc- 
tuations of thermal or quantum mechanical origin assist 
directed transport while large fluctuations (correspond- 
ing to high temperatures or to large values of ft e ff respec- 
tively) blur the asymmetry of the attractor and thus kill 
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Figure 2: (Color online) Asymptotic current value Joo as a 
function of temperature T. Upper panel shows the case for 
T = 0.7, the middle one for F = 0.75, and the lower one 
for r = 0.9. Black circles stand for the classical values. 
The quantum cases correspond to h = 0.055 (blue triangles), 
h — 0.165 (green diamonds), and h = 0.494 (red squares). 
As initial conditions we have taken 10 6 random points (clas- 
sical) and a density operator with equal population at all the 
possible p eigenstates inside the phase space region given by 
(x e [0,tt), p e [-tt,tt]). 

the net current. 

It is interesting to note that in the case F = 0.75 the 
thermal diffusion associated with the temperature which 
gives the maximal value of the current (< £ 2 >= 0.12 for 
T = 0.25) is of the order of the quantum coarse-graining 
ft e ff corresponding to the strongest quantum current. For 
T = 0.7 the classical current attains its maximum value 
at T = 0.05 (< £ 2 >= 0.03), which corresponds to a 
value of ft e ff that we were not able to consider in our 
numerical calculations. 

The analogy between thermal noise and quantum 
coarse-graining can also be appreciated when looking at 
the asymptotic Poincare sections and Husimi distribu- 
tions (displayed in Fig. 3 for T = 0.75, ft e ff = 0.055). 
As expected, at zero temperature the quantum Husimi 
function reproduces well the main patterns of the classi- 
cal attractor but shows less fine structure (see the upper 
panels). If a small temperature is introduced the fine de- 
tails of the classical distribution are smoothed out and 
both distributions look more alike (see the lower pan- 
els corresponding to T = 0.1). On the other hand the 
quantum distributions at zero and finite temperatures are 
practically indistinguishable, indicating that the quan- 
tum coarse-graining is at least of the order of the thermal 
one for these values of T. 

We finally study J(t) as a function of t (i.e., the number 
of iterations of the map). Results for T = 0.75 are shown 
in Fig. 4, where different temperatures and ft e ff values 
have been considered. Besides the mentioned fact that 
the asymptotic value is reached very rapidly, we notice 




Figure 3: Phase space portraits for F = 0.75 at t = 40. Left 
panels correspond to the classical while the right ones to the 
quantum strange attractors. In the upper panels T = 0, while 
in the lower ones T — 0.1. As initial conditions we have taken 
10 6 random points (classical) and a density operator with 
equal population at all the possible p eigenstates inside the 
phase space region given by (x £ [0,tt), p £ [— 7t,7t]). 



that the transient behavior shows a very close quantum- 
to-classical correspondence. The classical current peak 
observed at t ~ 10 for low temperatures (T = 0.1) is also 
present in the quantum current when ft e ff = 0.055 and 
T = 0. This peak disappears from the classical current at 
larger temperatures (T = 0.85), and so does the quantum 
one at ft e ff = 0.494 and T = 0. So the analogy seems to 
hold at all times. 



IV. CONCLUSIONS 

In this work we have analyzed the behavior of a chaotic 
dissipative system that shows directed transport under 
the influence of a thermal bath, both in its classical and 
quantum versions. We have varied the strength of the 
coupling with the environment and also the tempera- 
ture. We have found that the transport enhancement 
effect due to a finite temperature is highly dependent on 
the system specific properties. In fact, it depends on 
the coupling strength of the system with the environ- 
ment and also on the h size. Moreover, we could find 
an analogy between the effects caused by thermal and 
quantum fluctuations. These results open the possibility 
for many further studies that include finding ways of ob- 
taining large ratchet currents in experimentally realistic 
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situations in kicked BECs and cold atoms experiments. 
These are one of the best candidates to test our results 
since even BECs show an unavoidable fraction of non- 
condensed atoms when kicked. If kicks become strong, 
thermal excitations will be of much relevance rather than 
a negligible effect. With the sort of calculations presented 
in this paper the effects of this fraction on the transport 
properties of the system could be estimated. 



Figure 4: (Color online) Current J as a function of time t, 
for coupling strength T = 0.75. Thin black lines correspond 
to the classical values, while thick green (gray) ones to the 
quantum cases (in these last cases we always take T = 0) . We 
show results for T = 0.1 and ft e ff = 0.055 (dot-dashed lines), 
and for T = 0.85 and h e ft = 0.494 (solid lines). As initial 
conditions we have taken 10 7 random points (classical) and a 
density operator with equal population at all the possible p 
eigenstates inside the phase space region given by (x £ [0,7r), 
p e f-7r,7rl). 
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